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1 Introduction 

According to general relativity, not only the gravitational field but also the 
structure of space and time, the stage for all the other fields, is governed by 
the dynamical laws of physics. The space we see is not a fixed background, but 
it evolves on large time scales, even to such extreme situations as singularities 
where all of space collapses into a single point. At such a point, however, 
energy densities and tidal forces diverge; all classical theories break down, 
even general relativity itself. This implies that general relativity cannot be 
complete since it predicts its own breakdown. Already for a long time, it has 
been widely expected that a quantum theory of general relativity would cure 
this problem, providing a theory which can tell us about the fate of a classical 
singularity. 

Most of the time, quantum gravity has been regarded as being far away 
from observational tests. In such a situation, different approaches would have 
to be judged purely on grounds of internal consistency and their ability to 
solve conceptual problems. Those requirements are already very restrictive for 
the quantization of a complicated theory as general relativity, to the extent 
that in all the decades of intense research not a single completely convincing 
quantum theory of gravity has emerged yet, even though there is a number 
of promising candidates with different strengths. Still, the ultimate test of a 
physical theory must come from a confrontation with observations of the real 
world. For quantum gravity, this means observations of effects which happen 
at the smallest scales, the size of the Planck length £p as 10 _32 cm. 

In particular in the light of recent improvements in precision cosmology, 
the cosmological arena seems to be most promising for experimental tests. 
This is fortunate since also many conceptual issues arise in the cosmological 
setting where the universe is studied as a whole. Examples are the singularity 
problem mentioned above and the so-called problem of time which we will 
address later. Therefore, one can use the same methods and approximations 
to deal with conceptual problems and to derive observational consequences. 
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One of the main approaches to quantum gravity is based on a canonical 
quantization of general relativity, which started with the formal Wheeler- 
DeWitt quantization and more recently evolved into quantum geometry. Its 
main strength is its background independence, i.e. the metric tensor which 
describes the geometry of space is quantized as a whole and not split into a 
background and a dynamical part. Since most familiar techniques of quantum 
field theory rely on the presence of a background, for this ambitious approach 
new techniques had to be invented which are often mathematically involved. 
By now, most of the necessary methods have been developed and we are 
ready to explore them in simple but physically interesting situations. 

Being the quantization of a complicated, non-linear field theory, quantum 
gravity cannot be expected to be easily understood in full generality. As 
always in physics, one has to employ approximation techniques which isolate 
a small number of objects one is interested in without taking into account 
all possible interactions. Prominent examples are symmetric models (which 
are usually called minisuperspaces in the context of general relativity) and 
perturbations of some degrees of freedom around a simple solution. This 
opens the possibility to study the universe as a whole (which is homogeneous 
and isotropic at large scales) as well as the propagation of a single particle 
in otherwise empty space (where complicated interactions can be ignored). 

These two scenarios constitute the two main parts of this article. In the 
context of the first one (Section 5) we discuss the basic equations which 
govern the quantum evolution of an isotropic universe and special proper- 
ties which reflect general issues of quantum gravity. We then analyze these 
equations and see how the effects of quantum geometry solve and elucidate 
important conceptual problems. Quantum gravity effects in these regimes 
also lead to modifications of classical equations of motion which can be used 
in a phenomenological analysis. A different kind of phenomenology, related to 
the propagation of particles in empty space, is discussed in Section 6. In this 
context cosmological scales are involved for many proposals of observations, 
and so they fit into the present scheme. 

Both settings are now at a stage where characteristic effects have been 
identified and separated from the complicated, often intimidating technical 
foundation. This is a natural starting point for phenomenological analyzes 
and opens a convenient port of entry for beginners to the field. 

The article is intended to describe the basic formalism to an extent which 
makes it possible to understand the applications without requiring too much 
background knowledge (the presentation cannot be entirely background inde- 
pendent, though). The general framework of quantum geometry is reviewed 
briefly in Section 4 after recalling facts about general relativity (Section 2) 
and the Wheeler-DeWitt quantization (Section 3). For the details we provide 
a guide to the literature including technical reviews and original papers. 
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2 General relativity 

General relativity is a field theory for the metric g^ v on a space-time M which 
determines the line element 4 ds 2 = g^ lv (x)dx fl dx v . The line element, in turn 
specifies the geometry of space-time; we can, e.g., measure the length of a 
curve C:R-> M,t^ x»(t) using 

1(C) = Jds = J ^ 9llv {x{t))x»{t)x»{t)&t. 



2.1 Field equations 

While a space-time can be equipped with many different metrics, resulting 
in different geometries, only a subclass is selected by Einstein's field equa- 
tions of general relativity which are complicated non-linear partial differential 
equations with the energy density of matter as a source. 

They can be understood as giving the dynamical evolution of a space-like 
geometry in a physical universe. Due to the four-dimensional covariance, 
however, in general there is no distinguished space-like slice which could 
be used to describe the evolution. All possible slices are allowed, and they 
describe the same four-dimensional picture thanks to symmetries of the field 
equations. 

Selecting a slicing into space- like manifolds, the field equations take on 
different forms and do not show the four-dimensional covariance explicitly. 
However, such a formulation has the advantage that it allows a canonical 
formulation where the metric q ab only of space-like slices plays the role of co- 
ordinates of a phase space, whose momenta are related to the time derivative 
of the metric, or the extrinsic curvature K ab = — \q a b of a slice [1]. This is in 
particular helpful for a quantization since canonical quantization techniques 
become available. The momentum conjugate to the metric q ab is related to 
the extrinsic curvature by 

TT ab = -±^drtq~(K ab - q ab K c c ) 

where indices are raised by using the inverse q ab of the metric. The dynamical 
field equation, the analog of Einstein's field equations, takes the form of a 
constraint, 5 the Hamiltonian constraint 

8nG^detq(q ac q bd + q ad q bc - qabqcd)n ab ir cd - ^_ V dct 1 

+ yfdetq /9 mat tor(<?) = (1) 

4 In expressions with repeated indices, a summation over the allowed range is 
understood unless specified otherwise. We use greek letters fj,, v, . . . for space- 
time indices ranging from zero to three and latin letters from the beginning of 
the alphabet, a,b, . . . for space indices ranging from one to three. 

5 Note that this requires a relation between the basic fields in every point of space; 
there are infinitely many degrees of freedom and infinitely many constraints. 



4 



Martin Bojowald and Hugo A. Morales- Tecotl 



where G is the gravitational constant, 3 R(q) the so-called Ricci scalar of 
the spatial geometry (which is a function of the metric q), and /9 ma tter(<7) is 
the energy density of matter depending on the particular matter content (it 
depends on the metric, but not on its momenta in the absence of curvature 
couplings). 

The complicated constraint can be simplified slightly by transforming to 
new variables [2], which has the additional advantage of bringing general 
relativity into the form of a gauge theory, allowing even more powerful math- 
ematical techniques. In this reformulation, the canonical degrees of freedom 
are a densitized triad Ef which can be thought of as giving three vectors la- 
belled by the index 1 < i < 3. Requiring that these vectors are orthonormal 
defines a metric given by 

q ab = ^tEj\(E- 1 ) i a (E-y b . 

Its canonical conjugate is the Ashtekar connection 

K = K iK (2) 

where is the spin connection (given uniquely by the triad such that d a Ef + 
e ijkElE\ = 0) and K l a is the extrinsic curvature. The positive Barbero- 
Immirzi parameter 7 also appears in the symplectic structure together with 
the gravitational constant G 

{Ai(x),E b j (y)}=8n 7 G5 b J i j S(x,y) (3) 

and labels equivalent classical formulations. Thus, it can be chosen arbitrarily, 
but the freedom will be important later for the quantum theory. The basic 
variables can be thought of as a "vector potential" A\ and the "electric field" 
Ef of a gauge theory, whose gauge group is the rotation group SO (3) which 
rotates the three triad vectors: Ef \— > AjEj for A e SO (3) (such a rotation 
does not change the metric). 

Now, the Hamiltonian constraint takes the form [3]: 

I det Ef I -^djkF^EJEl - 2 ( 1 + 7 2 ) | dot £f | -^K^K^ EfE{ 

+8nG^\teUSf\ p mattC r(E) = (4) 

where F* b the curvature of the Ashtekar connection, and the matter energy 
density Pmatter (E) now depends on the triad via the metric. 

2.2 Approximations 

Given the complicated nature of the field equations, one has to resort to 
approximation schemes in order to study realistic situations. In the case of 
gravity, the most widely used approximations are: 
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Assume symmetries. This simplifies the field equations by eliminating 
several degrees of freedom and simplifying the relations between the re- 
maining ones. In a cosmological situation, for instance, one can assume 
space to be homogeneous such that the field equations reduce to ordinary 
differential equations in time. 

Perturbations around a simple known solution. One can, e.g., study a 
small amount of matter, e.g. a gravitational wave or a single particle, and 
its propagation in Minkowski space. To leading order, the back reaction of 
the geometry, which changes due to the presence of the particle's energy 
density, on the particle's propagation can be ignored. 
Asymptotic regimes with boundary conditions. In many situations it is 
possible to isolate interesting degrees of freedom by looking at boundaries 
of space-time with special boundary conditions capturing the physical 
situation. It can then be possible to ignore interactions with the bulk 
degrees of freedom which simplifies the analysis. This strategy is most 
widely used in the context of black hole physics, in its most advanced 
form with isolated horizon conditions; see, e.g. [4]. 

The first two approximation schemes and their applications in quantum 
geometry will be discussed on Sections 5 and 6, respectively. Since the last one 
so far does not have many cosmological applications, it will not be used here. 
It does have applications in quantum geometry, however, in the calculation 
of black hole entropy [5] . In this section we only illustrate the first one in the 
context of isotropic cosmology. 

2.3 Cosmology 

In the simplest case of a cosmological model we can assume space to be 
isotropic (looking the same in all its points and in all directions) which implies 
that one can choose coordinates in which the line element takes the form 



with the scale factor a(t) (the evolving "radius" of the universe). The con- 
stant k can take the values k = for a spatially flat model (planar), k = 1 for 
a model with positive spatial curvature (spherical), and k = — 1 for a model 
with negative spatial curvature (hyperbolic). Einstein's field equations re- 
strict the possible behavior of a(t) in the form of the Fricdmann equation 



Since also the matter density p(a) enters, we can find a(t) only if we specify 
the matter content. Common choices are "dust" with p(a) oc a~ 3 or "radia- 
tion" with p(a) oc a~ 4 (due to an additional red-shift factor), which describe 



d.s 2 = -di 2 + a(i) 2 ((l - kr 2 )dr 2 + r 2 (dd 2 + sintfdv? 2 )) 



(5) 



[6] 




(6) 
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the matter degrees of freedom collectively. After choosing the matter con- 
tent, we just need to solve an ordinary differential equation. For radiation in 
a spatially flat universe, e.g., all solutions are given by a(t) oc y/t — to where 
to is an integration constant. 

In a more complicated but also more fundamental way one can describe 
the matter by using additional matter fields 6 which enter via their Hamilto- 
nian (or total energy) . This results in a system of coupled ordinary differential 
equations, one for the scale factor and others for the matter fields. A common 
example in cosmology is a scalar <j> which has Hamiltonian 

H^a) = \a-hl + o?W{<j ) ) (7) 

with its potential W and the scalar momentum p^ = a 3 tp. Note that it is 
important to keep track of the a-dependence in cosmology since a is evolving; 
in the usual formulas for Hamiltonians on Minkowski space a does not appear. 

The Friedmann equation is now given by (6) with energy density p(a) = 
H c f > (a)/a 3 . Now, the right hand side depends explicitly on (p and p,p which 
both depend on time. Their evolution is given by the Hamiltonian equations 
of motion 

<fi = {<j>,H 4> }=p 4> /a 3 (8) 
p 4 , = {p^H 4 ,} = -a 3 W'{<l>) . (9) 

By using the first equation one can transform the second one into a second 
order equation of motion for (f>: 

$ = -3aa- 1 4>-W'(4>) (10) 

which in addition to the usual force term from the potential has a friction 
term proportional to the first derivative of <f>. The friction is strongest for a 
rapid expansion. 

When we come close to a = 0, the kinetic term usually dominates and 
even diverges when a = 0. This is problematic and leads to the singularity 
problem discussed in the following subsection. However, the divergence occurs 
only when p^ ^ for small a, so one could try to arrange the evolution of 
the scalar such that the divergence is avoided. In addition to suppressing 
the diverging kinetic term, we have the additional welcome fact that p^, ~ 
implies 4> ~ 0o — const. The right hand side of the Friedmann equation then 
becomes constant, (a/a) 2 w A = (16irG/3)W((f>o) for k = 0. Its solutions 
are given by a(t) oc exp(\/yTi) which describes an accelerated expansion, or 
inflation. Though motivated in a different way here, inflation is deemed to 
be an important ingredient in cosmological model building, in particular for 
structure formation. 

6 In a homogeneous model, matter "fields" are also described by a finite number 
of parameters only, e.g. a single one for a scalar <j>. 
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Unfortunately, however, it is very difficult to arrange the evolution of the 
scalar in the way described here; for - in addition to introducing a new field, 
the inflaton <f> - it requires very special scalar potentials and also initial values 
of the scalar. A common choice is a quadratic potential W{<p) — ^mcf) 2 (e.g., 
for chaotic inflation) which requires a very small m (a very flat potential) for 
inflation to take place long enough, and also a huge initial value 4>o pushing it 
up to Planck values. There is a plethora of models with intricate potentials, 
all requiring very special choices. 

Inflation in general is the term for accelerated expansion [7] , i.e. a > 0. It is 
not necessarily of the exponential form as above, but can be parameterized by 
different ranges of the so-called equation of state parameter w which needs to 
be less than w < — 3 for inflation. It can be introduced by a phenomenological 
a-dependence of the energy density, 

p(a) oc a- 3 ^ +1 ) . (11) 

Note, however, that this is in general possible only with a-dependent w except 
for special cases. Solutions for a (with k = 0) are then of the form 

{{t - i ) 2 /(3™+3) f or _ 1 < w < _ 1 (power-law inflation) 
exp(v^li) for w = — 1 (standard inflation) (12) 

{t - t) 2 /( 3 ™+ 3 ) for w < -1 (super-inflation) 

where t is an initial value (replaced by A — (16ir/3)Gp with the constant en- 
ergy density p for standard inflation). Note in particular that super-inflation 
(also called pole-law inflation) can be valid only during a limited period of 
time since otherwise a would diverge for t = to- While these possibilities add 
more choices for model building, they share with standard inflation that they 
arc difficult to arrange with scalar potentials. 

2.4 Singularities 

Trying to suppress the kinetic term has led us to introduce inflation as an 
ingredient in cosmological models. Can it lead to a regular evolution, provided 
we manage to arrange it in some way? The answer is no, for the following 
intuitive reason: We can get to be very small by making special choices, 
but it will not be exactly zero and eventually the diverging a~ 3 will win 
if we only go close enough to a = 0. In the end, we always have to face 
the singularity problem illustrated by the simple solution a(t) oc \/t — to for 
radiation: a(to) = such that all of space collapses to a single point (any 
length of a space-like curve at to measured with the line element (5) is zero) 
and the energy density diverges. The most dooming consequence is that the 
evolution breaks down: We cannot set up an initial value problem at to and 
evolve to values of t smaller than to- The theory does not tell us what happens 
beyond to- This consequence is a general property of general relativity which 
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cannot be avoided. We used the symmetric situation only for purposes of 
illustration, but the singularity problem remains true for any solution [8]. 
There will always be points which can be reached in a finite amount of time, 
but we will not be able to know anything as to what happens beyond such 
a point. General relativity cannot be complete since it predicts situations 
where it breaks down. 

This is the classical situation. Can it be better in a quantum theory of 
gravity? In fact, this has been the hope for decades, justified by the following 
motivation: The classical hydrogen atom is unstable, but we know well that 
quantum mechanics leads to a ground state of finite energy 
which cures the instability problem. One can easily see that this is the only 
non-relativistic energy scale which can be built from the fundamental param- 
eters purely for dimensional reasons. In particular, a non-zero h is necessary, 
for in the classical limit h — > the ground state energy diverges leading to the 
classical instability. As an additional consequence we know that the existence 
of a non-zero H leads to discrete energies. 

In gravity the situation is similar. We have its fundamental parameter G 
from which we can build a natural length scale 7 £p = V8ttGH, the Planck 
length. It is very tiny and becomes important only at small scales, e.g. close 
to classical singularities where the whole space is small. Where the Planck 
length becomes important we expect deviations from the classical behavior 
which will hopefully cure the singularity problem. In the classical limit, h — > 
0, the Planck length becomes zero and we would get back the singularity. 
Completing our suggestions from the hydrogen atom, we also expect discrete 
lengths in a quantum theory of gravity, the explicit form of which can only 
be concluded from a precise implementation. 

3 Wheeler-DeWitt quantum gravity 

As discussed, one can bring general relativity into a canonical formulation 
where the metric q a b and its momenta ir ab play the role of phase space coor- 
dinates (infinitely many because they depend on the points of space) , together 
with possible matter degrees of freedom and their momenta. This allows us 
to perform a canonical quantization (see, e.g., [9]) by representing quantum 
states as functionals ^(qab, 4>) of the metric and matter fields, corresponding 
to a metric representation. The metric itself then acts as a multiplication 
operator, and its conjugate ir ab by a functional derivative 7r ab = —ihd/dq a b- 
These arc the basic operators from which more complicated ones can be con- 
structed. 

7 Sometimes the Planck length is defined as VGH. 
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3.1 The Wheeler DeWitt equation 

In a canonical formulation of general relativity, the dynamics is determined 
by a constraint equation, (1) in the variables used here. Replacing q a b and TT ab 
by the respective operators yields a complicated constraint operator Hadm 
acting on a wave function $\ Since the classical expression must vanish, only 
states \P are allowed which are annihilated by the constraint operator, i.e. 
they have to fulfill the Whecler-DcWitt equation Hadm^ = 0. Since the 
constraint is quadratic in the momenta, this is a second order functional 
differential equation. However, it is only formal since it contains products 
of functional derivatives which have to be regularized in a way which does 
not spoil the properties of the theory, in particular its background indepen- 
dence. Such a regularization is complicated because the classical constraint 
is not even a polynomial in the basic fields, and so far it has not been done 
successfully in the ADM formulation. 

There is another apparent difficulty with the constraint equation: It is sup- 
posed to give us the dynamics, but there is no time dependence at all, and no 
time derivative part as in a Schrodinger equation. This is a general property 
of theories as general relativity which are invariant under four-dimensional 
coordinate transformations. We do not have an absolute notion of time, and 
thus it cannot appear in the basic evolution equation. Classically, we can 
introduce a time parameter (coordinate time t), but it just serves to pa- 
rameterize classical trajectories. It can be changed freely by a coordinate 
transformation. In the quantum theory, which is formulated in a coordinate 
independent way, coordinate time cannot appear explicitly. Instead, one has 
to understand the evolution in a relational way: there is no evolution with 
respect to an absolute time, but only evolution of all the degrees of freedom 
with respect to each other. After all, this is how we perceive time. We build 
a clock, which is a collection of matter degrees of freedom with very special 
interactions with each other, and observe the evolution of other objects, de- 
grees of freedom with weak interactions with the clock, with respect to its 
progression. Similarly, we can imagine to select a particular combination of 
matter or metric degrees of freedom as our clock variable and re-express the 
constraint equation as an evolution equation with respect to it [10, 11]. For 
instance, in a cosmological context we can choose the volume of space as 
internal time and measure the evolution of matter degrees of freedom with 
respect to the expansion or contraction of the universe. In general, however, 
a global choice of a time degree of freedom which would allow us to bring the 
full Wheeler-DeWitt equation into the form of an evolution equation, is not 
known; this is the problem of time in general relativity. 

Due to the complicated regularization and interprctational issues, appli- 
cations of the full Wheeler-DeWitt equation have been done only at a formal 
level for semiclassical calculations. 



10 



Martin Bojowald and Hugo A. Morales- Tecotl 



3.2 Minisuperspaces 

In order to study the theory explicitly, we again have to resort to approxi- 
mations. A common simplification of the Wheeler-DeWitt formalism is the 
reduction to minisupcrspace models where the space is homogeneous or even 
isotropic. Therefore, the metric of space is specified by a finite number of 
parameters only only the scale factor a in the isotropic case. While this 
is similar in spirit to looking for symmetric classical solutions as we did in 
section 2, there is also an important difference: If we want the symmetry to 
be preserved in time we need to restrict the time derivative of the metric, i.e. 
its canonical conjugate, in the same symmetric form. This is possible classi- 
cally, but in quantum mechanics it violates Heisenberg's uncertainty relations 
for the excluded degrees of freedom. Minisuperspace models do not just give 
us particular, if very special exact solutions as in the classical theory; their 
results must be regarded as approximations which are valid only under the 
assumption that the interaction with the excluded parameters is negligible. 

An isotropic minisupcrspace model has the two gravitational parameters 
a and its conjugate p a = 3aa/8irG together with possible matter degrees of 
freedom which we simply denote as <j) and p^. Using a Schrodinger quan- 
tization of the momenta acting on a wave function ip(a,<f>), the Fricdmann 
equation (6) is quantized to the Wheeler-DeWitt equation 



with matter Hamiltonian H^a). This equation is not unique due to ordering 
ambiguities on the left hand side. Here, we use the one which is related 
to the quantization derived later. Without fixing the ordering ambiguity, 
consequences derived from the equation are ambiguous [12]. 

The Wheeler-DeWitt equation quantizes the dynamical classical equation 
and thus should describe the quantum dynamics. As described before, in an 
isotropic model we can select the scale factor a as an internal time; evolution 
of the matter fields will then be measured not in absolute terms but in relation 
to the expansion or contraction of the universe. Interpreting a as a time 
variable immediately brings Eq. (13) to the form of a time evolution equation, 
albeit with an unconventional time derivative term. 

An unrcsolvablc problem of the Wheeler-DeWitt quantization, however, is 
that it is still singular. Energy densities, all depending on the multiplication 
operator a -1 are still unbounded, and the Wheeler-DeWitt equation does 
not tell us what happens at the other side of the classical singularity at 
a = 0. Instead, the point of view has been that the universe is "created" 
at a = such that initial conditions have to be imposed there. DeWitt 
[10] tried to combine both problems by requiring tp(0) = which can be 
interpreted as requiring a vanishing probability density to find the universe 
at the singularity. However, this very probability interpretation, which is 




(13) 
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just taken over from quantum mechanics, is not known to make sense in a 
quantum cosmological context. Furthermore, at the very least one would also 
need appropriate fall-off conditions for the wave function since otherwise we 
can still get arbitrarily close to the singularity. Appropriate conditions are 
not known, and it is not at all clear if they could always be implemented. The 
worst problem is, however, that DeWitt's initial condition is not well-posed 
in more general models where its only solution would vanish identically. 

DeWitt's condition has been replaced by several proposals which are mo- 
tivated from different intuitions [13, 14]. However, they do not eliminate the 
singularity; with them the wave function would not even vanish at a = in 
the isotropic case. They accept the singularity as a point of creation. 

Thus, the hope motivated from the hydrogen atom has not materialized. 
The isotropic universe model is still singular after quantizing. Do we have 
to accept that singularities of gravitational systems cannot be removed, not 
even by quantization? If the answer would be affirmative, it would spell severe 
problems for any desire to describe the real world by a physical theory. It 
would mean that there can be no complete description at all; any attempt 
would stop at the singularity problem. 

Fortunately, the answer turns out not to be affirmative. We will see in 
the following sections that singularities are removed by an appropriate quan- 
tization. Why, then, is this not the case for the Whccler-DeWitt quantiza- 
tion? One has to keep in mind that there is no mathematically well-defined 
Whccler-DeWitt quantization of full general relativity, and systematic inves- 
tigations of even the formal equations are lacking. What one usually does 
instead is merely a quantum mechanical application in a simple model with 
only a few degrees of freedom. There is no full theory which one could use to 
see if all quantization steps would also be possible there. General relativity is 
a complicated theory and its quantization can be done, if at all, only in very 
special ways which have to respect complicated consistency conditions, e.g. 
in the form of commutation relations between basic operators. In a simple 
model, all these problems can be brushed over and consistency conditions 
are easily overlooked. One hint that this in fact happened in the Wheelcr- 
DeWitt quantization is the lacking discreteness of space. We expected that a 
non-zero Planck length in quantum gravity would lead to the discreteness of 
space. While we did see the Planck length in Eq. (13), there was no associ- 
ated discreteness: the scale factor operator, which is simply the multiplication 
operator a, still has continuous spectrum. 

After the discussion it should now be clear how one has to proceed in 
a more reliable way. We have to use as much of the full theory of quantum 
gravity as we know and be very careful to use only techniques in our sym- 
metric models which can also be implemented in the full theory In this way, 
we would respect all consistency conditions and obtain a faithful model of 
the full theory. Ideally, we would even start from the full theory and define 
symmetric models there at the level of states and operators. 
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By now, we have good candidates for a full theory of quantum gravity, 
and in the case of quantum geometry [15, 16, 17] also a procedure to de- 
fine symmetric models from it [18]. We will describe the main results in the 
following two sections. 

4 Quantum geometry 

While the Wheeler-DeWitt quantization is based on the ADM formulation 
whose basic variables are the metric and extrinsic curvature of a spatial slice, 
quantum geometry is a newer approach based on Ashtekar's variables. Quan- 
tization, in particular of a non-linear field theory, is a delicate step whose 
success can depend significantly on the formulation of a given classical the- 
ory. Classical theories can usually be formulated in many different, equivalent 
ways, all being related by canonical transformations. Not all of these transfor- 
mations, however, can be implemented unitarily at the quantum level which 
would be necessary for the quantum theories to be equivalent, too. For in- 
stance, when quantizing one has to choose a set of basic variables closed under 
taking Poisson brackets which are promoted unambiguously to operators in 
such a way that their Poisson brackets are mapped to commutation relations. 
There is no unique prescription to quantize other functions on phase space 
which are not just linear functions of the basic ones, giving rise to quantiza- 
tion ambiguities. In quantum mechanics one can give quite general conditions 
for a representation of at least the basic variables to be unique (this represen- 
tation is the well-known Schrodinger quantization). However, such a theorem 
is not available for a field theory with infinitely many degrees of freedom such 
that even the basic variables cannot be quantized uniquely without further 
conditions. 

One can often use symmetry requirements together with other natural 
conditions in order to select a unique representation of the basic variables, 
e.g. Poincare invariance for a field theory on Minkowski space as a back- 
ground [19]. For general relativity, which is background independent, it has 
recently been proven in the context of quantum geometry that diffeomor- 
phism invariance, i.e. invariance under arbitrary deformations of space, can 
replace Poincare invariance in strongly restricting the class of possible repre- 
sentations [20]. It is clear that those precise theorems can only be achieved 
within a theory which is mathematically well-defined. The Wheeler-DeWitt 
quantization, on the other hand, does not exist beyond a purely formal level 
and it is unknown if it can give a well-defined quantum representation of the 
ADM variables at all. In any case, it is based on basic variables different from 
the ones quantum geometry is based on so that any representation it defines 
would likely be inequivalent to the one of quantum geometry. 

From the beginning, quantum geometry was striving for a mathematically 
rigorous formulation. This has been possible because it uses Ashtekar's vari- 
ables which bring general relativity into the form of a gauge theory. While 
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not all standard techniques for quantizing a gauge theory can be applied 
(most of them are not background independent), new powerful techniques 
for a background independent quantization have been developed [21, 22, 23]. 
This was possible only because the space of connections, which is the config- 
uration space of quantum geometry, has a structure much better understood 
than the configuration space of the Wheeler-DeWitt quantization, namely 
the space of metrics. 

We do not describe those techniques here and instead refer the interested 
reader to the literature where by now several technical reviews are available 
[17, 24]. In this section, instead, we present an intuitive construction which 
illustrates all the main results. 

4.1 Basic operators and states 

As usually in gauge theories (for instance in lattice formulations), one can 
form holonomies as functions of connections for all curves e: [0, 1] — > £ in a 
manifold S, 

h e (A) = V exp ^1 Ai(e{t)) e a (t) n d?j e SU(2) (14) 

where e a is the tangent vector to the curve e and Tj = — ^<7j are generators 
of the gauge group SU(2) in terms of the Pauli matrices. The symbol V 
denotes path ordering which means that the non-commuting su(2) elements 
in the exponential arc ordered along the curve. Similarly, given a surface 
S: [0, l] 2 — > £ we can form a flux as a function of the triads, 

E(S) = f E?(y)n a ( y yd 2 y (15) 
Js 

where n a is the co-normal 8 to the surface S. Holonomies and fluxes are the 
basic variables which are used for quantum geometry, and they represent the 
phase space of general relativity faithfully in the sense that any two config- 
urations of general relativity can be distinguished by evaluating holonomies 
and fluxes in them. 

One can now prove that the set of holonomies and fluxes is closed under 
taking Poisson brackets and that there is a representation of this Poisson 
algebra as an operator algebra on a function space. Moreover, using the ac- 
tion of the diffeomorphism group on U, which deforms the edges and surfaces 
involved in the above definitions, this representation is the unique covariant 
one. Note that unlike the functional derivatives appearing in the Whceler- 
DeWitt quantization, these are well-defined operators on an infinite dimen- 
sional Hilbert space. Note in particular that holonomies are well-defined as 

8 The co-normal is defined as n a — ^e abc e de (dx b / dy d )(dx c / dy e ) without using a 
background metric, where x a are coordinates of E and y d coordinates of the 
surface S. 
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operators, but not the connection itself. A Wheeler-DeWitt quantization, on 
the other hand, regards the extrinsic curvature, related to the connection, as 
one of the basic fields and would try to promote it to an operator. This is 
not possible in quantum geometry (and it is not known if it is possible at a 
precise level at all) which demonstrates the inequivalence of both approaches. 
The fact that only the holonomies can be quantized can also be seen as one of 
the consistency conditions of a full theory of quantum gravity mentioned ear- 
lier. In a minisupcrspace model one can easily quantize the isotropic extrinsic 
curvature, which is proportional to p a /a. However, since it is not possible in 
the full theory, the model departs from it already at a very basic level. A 
reliable model of a quantum theory of gravity should implement the feature 
that only holonomies can be quantized; we will come back to this issue later. 

We did not yet specify the space of functions on which the basic opera- 
tors act in the representation of quantum geometry. Understandably, a full 
definition involves many techniques of functional analysis, but it can also be 
described in intuitive terms. As mentioned already, it is convenient to define 
the theory in a connection representation since the space of connections is 
well-understood. We can then start with the function 1 which takes the value 
one in every connection and regard it as our ground state. 9 The holonomies 
depend only on the connection and thus act as multiplication operators in 
a connection formulation [25]. Acting with a single holonomy h e (A) on the 
state 1 results in a state which depends on the connection in a non-trivial 
way, but only on its values along the curve e. More precisely, since holonomies 
take values in the group SU(2), we should choose an SU(2)-representation, for 
instance the fundamental one, and regard the matrix elements of the holon- 
omy in this representation as multiplication operators. This can be done with 
holonomies along all possible curves, and also acting with the same curve sev- 
eral times. Those operators can be regarded as basic creation operators of 
the quantum theory. Acting with holonomies along different curves results 
in a dependence on the connection along all the curves, while acting with 
holonomies along the same curve leads to a dependence along the curve in 
a more complicated way given by multiplying all the fundamental represen- 
tations to higher ones. One can imagine that the state space obtained in 
this way with all possible edges (possibly intersecting and overlapping) in 
arbitrary numbers is quite complicated, but not all states obtained in this 
way are independent: one has to respect the decomposition rules of repre- 
sentations. This can all be done resulting in a basis of states, the so-called 
spin network states [26]. Furthermore, they are orthonormal with respect to 

9 Note that we do not call it "vacuum state" since the usual term "vacuum" 
denotes a state in which matter is unexcited but the gravitational background is 
Minkowski space (or another non-degenerate solution of general relativity). We 
will see shortly, however, that the ground state we are using here represents a 
state in which even gravity is "unexcited" in the sense that it defines a completely 
degenerate geometry. 
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the diffcomorphism invariant measure singled out by the representation, the 
Ashtekar-Lewandowski measure [22] . 

Note also that the quantum theory should be invariant under SU(2)- 
rotations of the fields since a rotated triad does not give us a new metric. 
Holonomies are not gauge invariant in this sense, but as in lattice gauge theo- 
ries we can use Wilson loops instead which are defined as traces of holonomies 
along a closed loop. Repeating the above construction only using Wilson loops 
results in gauge invariant states. 

Since we used holonomies to construct our state space, their action can be 
obtained by multiplication and subsequent decomposition in the independent 
states. Fluxes, on the other hand, are built from the conjugate of connections 
and thus become derivative operators. Their action is most easy to under- 
stand for a flux with a surface S which is transversal to all curves used in 
constructing a given state. Since the value of a triad in a given point is conju- 
gate to the connection in the same point but Poisson commutes with values 
of the connection in any other point, the flux operator will only notice inter- 
section points of the surface with all the edges which will be summed over 
with individual contributions. The contributions of all the intersection points 
are the same if we count intersections with overlapping curves separately. In 
this way, acting with a flux operator on a state returns the same state multi- 
plied with the intersection number between the surface of the flux and all the 
curves in the state. This immediately shows us the eigenvalues of flux oper- 
ators which turn out to be discrete. Since the fluxes are the basic operators 
representing the triad from which geometric quantities like length, area and 
volume are constructed, it shows that geometry is discrete [27, 28, 29, 30]. The 
main part of the area spectrum for a given surface S (the one disregarding 
intersections of curves in the state) is 

i4(S)= 7 4£y*(ji + l) (16) 

i 

where the sum is over all intersections of the surface S with curves in the 
state, and the SU(2)-labcls ji parameterize the multiplicity if curves overlap 
(without overlapping curves, all ji are |). Thus, quantum geometry predicts 
that geometric spectra are discrete, and it also provides an explicit form. 
Note that the Planck length appears (which arises because the basic Poisson 
brackets (3) contain the gravitational constant while h enters by quantizing a 
derivative operator), but the scale of the discreteness is set by the Barbero- 
Immirzi parameter 7. While different 7 lead to equivalent classical theories, 
the value of the parameter does matter in the quantum theory. If 7 would be 
large the discreteness would be important already at large scales despite the 
smallness of the Planck length. Calculations from black hole entropy, however, 
show that 7 must be smaller than one, its precise value being log(2)/7T\/3 [5]. 

Thus, quantum geometry already fulfills one of our expectations of Sec- 
tion 2, namely that quantum gravity should predict a discreteness of geometry 
with a scale set roughly by the Planck length. Note that the use of holonomies 
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in constructing the quantum theory, which was necessary for a well-defined 
formulation, is essential in obtaining the result about the discreteness. This 
fact has been overlooked in the Wheeler-DeWitt quantization which, conse- 
quently, does not show the discreteness. 



4.2 Composite operators 

We now have a well-defined framework with a quantization of our basic quan- 
tities, holonomies and fluxes. Using them we can construct composite oper- 
ators, e.g. geometric ones like area and volume or the constraint operator 
which governs the dynamics. Many of them have already been defined, but 
they are quite complicated. The volume operator, for instance, has been con- 
structed and it has been shown to have a discrete spectrum [27, 29, 31]; 
determining all its eigenvalues, however, would require the diagonalization 
of arbitrarily large matrices. Since it plays an important role in constructing 
other operators, in particular the Hamiltonian constraint [32, 33], it makes 
explicit calculations in the whole theory complicated. 

The constraint, for instance can be quantized by using a small Wilson 
loop along some loop a, which has the expansion h a = 1 + As^s^F^Ti where 
A is the coordinate area of the loop and si and S2 are tangent vectors to 
two of its edges, to quantize the curvature of the connection. The product 
of triads divided by the determinant appears to be problematic because the 
triad can be degenerate resulting in a vanishing determinant. However, one 
can make use of the classical identity [32] 

E a E*?e ijk 1 

' 1 = -r^^ abc {A k c ,V} , (17) 
^det^ 47r 7 G c ' 

replace the connection components by holonomies h s and use the volume 
operator to quantize this expression in a non-degenerate way. For the first 
part 10 of the constraint (4) this results in 

ff=E E e IJK tv(h aij{A) h SK{A) [h^ {A) ,V v ]) (18) 

vEV v(A)=v 

where we sum over the set V of vertices of the graph belonging to the state 
we act on, and over all possible choices (up to diffeomorphisms) to form a 
tetrahedron A with a loop au(A) sharing two sides with the graph and a 
third transversal curve Sk{A). The first holonomy along au(A) quantizes the 
curvature components while h SK (a) together with the commutator quantizes 
the triad components. 



10 The remaining part of the constraint involving extrinsic curvature components 
can be obtained from the first part since the extrinsic curvature can be written 
as a Poisson bracket of the first part of the constraint with the volume [32]. 
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A similar strategy can be used for matter Hamiltonians which usually 
also require to divide by the determinant of the triad at least in their kinetic 
terms. Here it is enough to cite some operators as examples to illustrate the 
general structure which will later be used in Section 6; for details we refer to 
[33]. For electromagnetism, for instance, we need to quantize 



-^Maxwell 



Qab 



,[E a E b + B a B b ] 



(19) 



2Q 2 Vdetg 

where Q4 a ,£: a /Q 2 ) are the canonical fields of the electromagnetic sector with 
gauge group U(l) and coupling constant Q, related to the dimensionless fine 
structure constant by ckem = Q 2 h. Furthermore, i? b is the magnetic field of 
the U(l) connection A, i.e. the dual of its curvature. 

Along the lines followed for the gravitational Hamiltonian we obtain the 
full electromagnetic Hamiltonian operator [33] (with a weight factor w(v) 
depending on the graph) 



H 



Maxwell 



E w ^ 



■>t '()-> " " ■■ E 

p ^ veV v(A)=v(A>) 



tr n h t 



l(A) 



x tr Tih Sp{AI) 



sr(A>)> 



w JKL MNP 



(e-i^K _ ^-^m- - l) - 



V JK V MN) 



(20) 

Let us emphasize the structure of the above regularized Hamiltonian. There 
is a common gravitational factor included in the SU(2) trace. The basic en- 
tities that quantize the electromagnetic part are the corresponding fluxes (as 
operators acting on a state for the electromagnetic field): one is associated 
with the magnetic field, which enters through a product of exponential flux 
factors exp(— i^')- 8 ) constructed from holonomies in A and A', respectively, 
while the other is related to the electric field, entering in a bilinear product 
of electric fluxes . 

Similarly, one can set up a theory for fermions which would be coupled to 
the gauge fields and to gravity. For a spin-i field 9 a one obtains the kinetic 
part (with the Planck mass mp = h/lp) 



H, 



spin— 1/2 



mp 

M 



E E 

veVv(A)=' 



x tr Tjh sj(A) 



e ijk e UK tr 



sj(A) ' 



Tih 



[{Tkh SK (A)Q)\s K (A) - 0\v]_ 



si (A) 



s,(A)> 



h.c. 



(21) 



lA d6 A (v) 

All of these matter Hamiltonians are well-defined, bounded operators, 
which is remarkable since in quantum field theories on a classical background 
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matter Hamiltonians usually have ultraviolet divergences. This can be inter- 
preted as a natural cut-off implied by the discrete structure. Compared to the 
Whceler-DeWitt quantization it is a huge progress that well-defined Hamil- 
tonian constraint operators are available in the full theory. Not surprisingly, 
their action is very complicated for several reasons. The most obvious ones 
are the fact that Wilson loops necessary to quantize curvature components 
create many new curves in a state which is acted on, and that the volume 
operator is being used to quantize triad components. The first fact implies 
that complicated graphs are created, while the second one shows that even a 
single one of those contributions is difficult to analyze due to the unknown 
volume spectrum. And after determining the action of the constraint opera- 
tor on states we still have to solve it, i.e. find its kernel. Furthermore, there 
are always several possible ways to quantize a classical Hamiltonian such 
that the ones we wrote down should be considered as possible choices which 
incorporate the main features. 

The complicated nature should not come as a surprise, though. After 
all, we are dealing with a quantization of full general relativity without any 
simplifying assumptions. Even the classical equations are difficult to solve and 
to analyze if we do not assume symmetries or employ approximation schemes. 
Those simplifications are also available for quantum geometry, which is the 
subject of the rest of this article. Symmetries can be introduced at the level of 
states which can be rigorously defined as distributional, i.e. non-normalizable 
states (they cannot be ordinary states since the discrete structure would break 
any continuous symmetry) . Approximations can be done in many ways, and 
different schemes are currently being worked out. 

5 Loop quantum cosmology 

Loop quantum cosmology aims to investigate quantum geometry in simplified 
situations which are obtained by implementing symmetries. In contrast to a 
Whceler-DeWitt quantization and its minisuperspace models, there is now 
also a full theory available. It is then possible to perform all the steps of the 
quantization in a manner analogous to those in the full theory. In particular, 
one can be careful enough to respect all consistency conditions as, e.g., the 
use of holonomies. 11 There is a tighter relation between symmetric models 
and the full theory which goes beyond pure analogy. For instance, symmetric 
states can be defined rigorously [18, 34, 35] and the relation between operators 
is currently being investigated. In this section, as already in the previous one, 
we use intuitive ideas to describe the results. 

11 This sometimes requires to perform manipulations which seem more complicated 
than necessary or even unnatural from the point of view of a reduced model. 
However, all of this can be motivated from the full theory, and in fact exploiting 
simplifications which are not available in a full theory can always be misleading. 
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In addition to testing implications of the full theory in a simpler context, 
it is also possible to derive physical results. Fortunately, many interesting and 
realistic physical situations can be approximated by symmetric ones. This is 
true in particular for cosmology where one can assume the universe to be 
homogeneous and isotropic at large scales. 



5.1 Symmetric states and basic operators 

As seen before, the canonical fields of a theory of gravity are completely 
described by two numbers (depending on time) in an isotropic context. For 
Ashtekar's variables, isotropic connections and triads take the form 

Ai(x)dx a = clo* , E?-^=pX i (22) 

where uf are invariant 1-forms and Xi invariant vector fields. For a spatially 
flat configuration, w ! = dx l are just coordinate differentials, while Xi are the 
derivatives; for non-zero spatial curvature the coordinate dependence is more 
complicated. In all isotropic models, c and p are functions just of time. Their 
relation to the isotropic variables used before is 

c=i(fc- 7 <x) , \p\=a 2 (23) 

such that {c,p} — (8ir/3)-fG. An important difference is that p can have both 
signs, corresponding to the two possible orientations of a triad. In a classical 
theory we would ultimately have to restrict to one sign since p = represents 
a degenerate triad, and positive and negative signs are disconnected. But the 
situation can (and will) be different in a quantum theory. 

We can now perform an analog of the construction of states in the full 
theory. The symmetry condition can be implemented by using only invariant 
connections (22) in holonomies as creation operators, i.e. 

hi(c) = cxp(cTi) = cos(c/2) + 2t{ sin(c/2) . 

Consequently, all the states we construct by acting on the ground state 1 
will be functions of only the variable c. All the complication of the full states 
with an arbitrary number of curves has collapsed because of our symmetry 
assumption. The analog of the spin network basis, an orthonormal basis in 
the connection representation, is given by 12 [37] 

= cxp(W2) 
V ' ' V2sin(c/2) V ' 

for all integer n. 

12 A more careful analysis shows that the Hilbert space of loop quantum cosmology 
is not separable [36]. For our purposes, however, we can restrict to the separable 
subspace used here which is left invariant by our operators. 
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An analog of the flux operator is given by a quantization of the isotropic 
triad component p, 

p\n) = \~flln\n) . (25) 

This immediately allows a number of observations: Its spectrum is discrete 
(n corresponds to the intersection number in the full theory) which results in 
a discrete geometry. The scale factor a = \/\p\ also has a discrete spectrum 
which is very different from the Wheeler-DeWitt quantization where the scale 
factor is just a multiplication operator with a continuous spectrum. Thus, we 
obtain a different quantization with deviations being most important for small 
scale factors, close to the classical singularity. The quantization p also tells 
us that the sign of the "intersection number" n determines the orientation of 
space. There is only one state, |0), which is annihilated by p; we identify it 
with the classical singularity. 

We can also use p in order to obtain a quantization of the volume [38]; 
here we use the convention 

V(H-i)/2 = (g74) 1 V(H-l)\n\(\n\ + l) (26) 

for its eigenvalue on a state |n). Thus, we achieved our aims: We have a 
quantization of a model which simplifies the full theory. This can be seen 
from the simple nature of the states and the explicit form of the volume 
spectrum which is not available in the full theory. As we will see later, this 
also allows us to derive explicit composite operators [39, 37]. At the same 
time, we managed to preserve essential features of the full theory leading 
to a quantization different from 13 the Wheeler-DeWitt one which lacks a 
relation to a full theory. 

5.2 Inverse powers of the scale factor 

We can now use the framework to perform further tests of essential aspects 
of quantum cosmology. In the Wheeler-DeWitt quantization we have seen 
that the inverse scale factor a -1 becomes an unbounded operator. Since its 
powers also appear in matter Hamiltonians, their quantizations will also be 
unbounded, reflecting the classical divergence. The situation in loop quan- 
tum cosmology looks even worse at first glance: a still contains zero in its 
spectrum, but now as a discrete point. Then its inverse does not even exist. 
There can now be two possibilities: It can be that we cannot get a quantiza- 
tion of the classically diverging a -1 , which would mean that there is no way 
to resolve the classical singularity. As the other possibility it can turn out 
that there are admissible quantizations of a~ l in the sense that they have 
the correct classical limit and are densely defined operators. The second pos- 
sibility exists because, as noted earlier, there are usually several possibilities 

13 Both quantizations are in fact inequivalent because, e.g., the operator a has 
discrete spectrum in one and a continuous spectrum in the other quantization. 
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to construct a non-basic operator like a -1 . If the simplest one fails (looking 
for an inverse of a), it does not mean that there is no quantization at all. 

It turns out that the second possibility is realized [40] , in a way special to 
quantum geometry. We can use the identity (17), which has been essential in 
quantizing Hamiltonians, in order to rewrite a -1 in a classically equivalent 
form which allows a well-defined quantization. In this way, again, we stay 
very close to the full theory, repeating only what can be done there, and at 
the same time obtain physically interesting results. 

The reformulation can be written in a simple way for a symmetric context, 



e.g., 14 



a = 



V2tt 7 G l V \37r7G 



= i^2i, 1 G ] { ] + l)(2 ] + l))- 1 ^ ] {rMc){h l {cr\^}^ (27) 

indicating in the last step that we can choose any SU(2)-representation when 
computing the trace without changing the classical expression (the trace with- 
out label is in the fundamental representation, j — |). Note that we only need 
a positive power of p at the right hand side which can easily be quantized. We 
just have to use holonomy operators and the volume operator, and turn the 
Poisson bracket into a commutator. This results in a well-defined operator 
which has eigenstates |n) and, for j = |, eigenvalues 

«^>> = (\AWs - v ^"l/ 2 - 1 ) 2 |n> (28) 



in terms of the volume eigenvalues (26). 

It has the following properties [40, 41]; see Fig. 1: 



It is a finite operator with upper bound 15 (a _1 ) max = 32 ^f^~ ^ at 
peak at n = 2. Now we can see that the situation is just as in the case 
of the hydrogen atom: The classically pathological behavior is cured by 
quantum effects which purely for dimensional reasons - require h in the 
denominator (recall £p = \/8ttGE). A finite value for the upper bound is 
possible only with non-zero £p; in the classical limit £p — > we reobtain 
the classical divergence. 



a 



One can easily see that there are many ways to rewrite a -1 in such a way. 
Essential features, however, are common to all these reformulations, for instance 
the fact that always the absolute value of p will appear in the Poisson bracket 
rather than p itself. We will come back to this issue shortly in the context of 
quantization ambiguities. 

We only use the general form of the upper bound. Its precise value depends on 
quantization ambiguities and is not important in this context. 
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n 



Fig. 1. Two examples for eigenvalues of inverse scale factor operators with j = 5 
(+) and j = 10 (x), compared to the classical behavior (dashed) and the approxi- 
mations (30) [41]. 

2. The first point demonstrates that the classical behavior is modified at 
small volume, but one can see that it is approached rapidly for volumes 
larger than the peak position. Thus, the quantization (28) has the correct 
classical limit and is perfectly admissible. 

3. While the first two points verify our optimistic expectations, there is also 
an unexpected feature. The classical divergence is not just cut off at a 
finite value, the eigenvalues of the inverse scale factor drop off when we go 
to smaller volume and are exactly zero for n = (where the eigenvalue of 
the scale factor is also zero). This feature, which will be important later, 
is explained by the fact that the right hand side of (27) also includes 
a factor of sgn(p) 2 since the absolute value of p appears in the Poisson 
bracket. Strictly speaking, we can only quantize sgn(p) 2 a _1 , not just a -1 
itself. Classically, we cannot distinguish between both expressions - both 
are equally ill-defined for a = and we would have to restrict to positive 
p. As it turns out, however, the expression with the sign does have well- 
defined quantizations, while the other one does not. Therefore, we have 
to use the sign when quantizing expressions involving inverse powers of a, 
and it is responsible for pushing the eigenvalue of the inverse scale factor 
at n = to zero. 

4. As already indicated in (27), we can rewrite the classical expression in 
many equivalent ways. Quantizations, however, will not necessarily be 
the same. In particular, using a higher representation j ^ i in (27), the 
holonomies in a quantization will change n by amounts larger than one. 
In (28) we will then have volume eigenvalues not just with n— 1 and n+1, 
but from n — 2j to n + 2j corresponding to the coupling rules of angular 
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momentum. Quantitative features depend on the particular value of j (or 
other quantization ambiguities), but qualitative aspects - in particular 
the ones in points 1 to 3 - do not change. Thus, the quantization is robust 
under ambiguities, but there can be small changes depending on which 
particular quantization is used. Such a freedom can also be exploited in 
a phenomenological analysis of some effects. 

Let us make the last point more explicit. The exact formula for eigenvalues 
with a non-fundamental representation is quite complicated. It can, however, 
be approximated using a rather simple function [41] 




Fig. 2. The function p(q) of (29), derived in [41]. For small q, p(q) increases like 
p(q) = fq r ' /4 (l-q + 0(q 2 )) (dashed). 



p(q) 



77 q 



7((g + l) 11/4 -k-i| 11/4 ) 

llg((g + l) 7 / 4 -sgn(g-l)|(,-l| 7 / 4 ) 



see Fig. 2, such that the eigenvalues of a quantization of a 
m are given by 

(«" m )i j) = ^i/fp(M/2j) 2m 



(29) 

with positive 
(30) 



with the ambiguity parameter j. There are many other ambiguities which can 
change also the function p, but the one indicated by j is most important. It 
parameterizes the position of the peak in an inverse power of the scale factor, 
which roughly coincides with the boundary between classical behavior and 
quantum modifications. 
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Note also that (30) displays the observation that inverse powers of the 
scale factor annihilate the singular state |0), thanks to ]xm q _>o(q~ 1 / 4, p(q)) = 
0. This has important consequences for matter Hamiltonians: they usually 
consist of a kinetic term containing components of the inverse metric and a 
potential term containing metric components. In the isotropic case, a widely 
used example in cosmology is the scalar Hamiltonian (7). We have already 
discussed the divergence at the singularity of the classical kinetic term, unless 
Ptf, — 0. The potential term, on the other hand, vanishes there. When we 
quantize this expression, we have to use an inverse power (30) in the kinetic 
term; for the potential term we can just use volume eigenvalues. Now, the 
potential term still vanishes at n = 0, but so does the kinetic term after 
quantization. Similarly, one can check that any matter Hamiltonian /f ma ttor 
(assuming for simplicity the absence of curvature couplings) fulfills 

ffmattcr|0) =0 (31) 

irrespective of the particular kind of matter and its quantization. 

All these observations indicate that the quantum behavior is much better, 
less singular, than the classical one. The real test, however, can only come 
from studying the quantum evolution. An absence of singularities can be 
confirmed only if it is possible to extend the evolution through the singular 
boundary; the theory has to tell us what happens at the singularity and 
beyond. 

5.3 Dynamics 16 

To study the dynamics of a theory we need its evolution equation which 
for gravity is given by the Hamiltonian constraint. In the Wheeler-DeWitt 
quantization we have seen that the constraint equation takes the form of an 
evolution equation after quantizing in a metric or triad representation and 
choosing an internal time a. 

We can follow the same steps here if we first transform from the con- 
nection representation used so far in quantum geometry to a triad repre- 
sentation. This can be done straightforwardly since we already know the 
triad eigenstates |n). A state \ip) can then be expanded in these eigenstates, 
\ip) — J2 n V'n(0)l n ) denoting possible matter degrees of freedom collectively 
by (j). The coefficients Tp n {4>) m the expansion then define, as usually, the 
state in the triad representation. Since n denotes the eigenvalues of p, it will 
now play the role of an internal time. Here we observe another difference to 
the Wheeler-DeWitt quantization: due to the discrete geometry, also time is 
discrete in an internal time picture. 

The Wheeler-DeWitt quantization now proceeded by quantizing the grav- 
itational momentum p a by a differential operator as in quantum mechanics. 

16 A brief summary of the results in this subsection can be found in [42]. 
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An analogous step is not possible in quantum geometry; momenta here, i.e. 
connection components, have to be quantized using holonomics which do not 
act as differential operators. Instead, they act according to the SU(2) coupling 
rules, e.g. 

(c\hi(c)\n) = (c|cos(c/2) + 2T i sin(c/2)|n) 

= ±«c|n + 1) + (c\n - 1)) - \in{{c\n + 1) - (c\n - 1)) . 

Thus, in a triad representation holonomies act by changing the label n in 

i>n{<t>) by ±1 since, e.g., 

(sin(c/2)V0„ = -\\^i) n {\n + 1) - \n - 1)) = ±i^(V> n +i - tp n -i)\n) . 

n n 

The constraint operator contains several holonomy operators and also the 
volume operator. It leads to the constraint equation [37, 43] 

(V\ n+4 \/2 - V|„ +4 |/2-l)c ife Vr l +4(</') - (2 + J 2 k 2 )(V lnl/2 - Vf n | /2 _l)^„(<£) 
+ (V|„-4|/2 - V| rl _4|/2-l)e _ife V'r l -4(0) 

= - T G7 3 4^matter(n)Vn(^) (32) 

which is a difference equation rather than a differential equation thanks to the 
discrete internal time. The parameter k again signifies the intrinsic curvature; 
for technical reasons the above equation has only been derived for the values 
k = and k = 1, not for k = —1. 

While the left hand side is very different from the Wheeler-DeWitt case, 
the right hand side looks similar. This is, however, only superficially so; for we 
have to use the quantizations of the preceding subsection for inverse metric 
components, in particular in the kinetic term. 

We can eliminate the phase factors e ±lfc in (32) by using a wave function 
V>7i(</>) := e mk ^^ n {4>) which satisfies the same equation without the phase 
factors (of course, it is different from the original wave function only for 
k = 1). The phase factor can be thought of as representing rapid oscillations 
of the wave function caused by non-zero intrinsic curvature. 

Large volume behavior 

Since the Wheeler-DeWitt equation corresponds to a straightforward quan- 
tization of the model, it should at least approximately be valid when we are 
far away from the singularity, i.e. when the volume is large enough. To check 
that it is indeed reproduced we assume large volume, i.e. n>l, and that the 
discrete wave function - 0n(</ ) ) does not display rapid oscillations at the Planck 
scale, i.e. from n to n+1, because this would indicate a significantly quantum 
behavior. We can thus interpolate the discrete wave function by a continuous 
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one tp(p, (f>) = V> n ( p )(0) with n(p) = 6p/j£p from (25). By our assumption of 
only mild oscillations, tp{p, <j>) can be assumed to be smooth with small higher 
order derivatives. We can then insert the smooth wave function in (32) and 
perform a Taylor expansion of tp n ±i{<t>) = ^{p( n ) ± §7^p) m terms of p/j£p- 
It is easy to check that this yields to leading order the equation 

\ (^^? ~ fc ) i>(P^) = -yGi? matt cr(p)V>(p,0) 

which with a = \/\p\ is the Wheeler-DeWitt equation (13) in the ordering 
given before [44]. Thus, indeed, at large volume the Wheeler-DeWitt equation 
is reproduced which demonstrates that the difference equation has the correct 
continuum limit at large volume. It also shows that the old Wheeler-DeWitt 
quantization, though not the fundamental evolution equation from the point 
of view of quantum geometry, can be used reliably as long as only situations 
are involved where the discreteness is not important. This includes many 
semiclassical situations, but not questions about the singularity. 

When the volume is small, we are not allowed to do the Taylor expansions 
since n is of the order of one. There we expect important deviations between 
the difference equation and the approximate differential equation. This is 
close to the classical singularity, where we want corrections to occur since 
the Wheeler-DeWitt quantization cannot deal with the singularity problem. 



Non-singular evolution 

To check the issue of the singularity in loop quantum cosmology we have to 
use the exact equation (32) without any approximations [45]. We start with 
initial values for tp n (<fi) at large, positive n where we know that the behavior is 
close to the classical one. Then, we can evolve backwards using the evolution 
equation as a recurrence relation for ipn-4,{<fi) m terms of the initial values. In 
this way, we evolve toward the classical singularity and we will be able to see 
what happens there. The evolution is unproblematic as long as the coefficient 
V\ n -4\/2 — V\n-4\/2-i of ip n -4 in the evolution equation is non-zero. It is easy 
to check, however, that it can be zero, if and only if n = 4. When n is four, 
we are about to determine the value of the wave function at n = 0, i.e. right 
at the classical singularity, which is thus impossible. It seems that we are 
running into a singularity problem again: the evolution equation does not 
tell us the value - 0o(0) there. 

A closer look confirms that there is no singularity. Let us first ignore the 
values ^0(0) and try to evolve through the classical singularity. First there 
are no problems: for tp_i(<p) we only need tp3(4>) and ipr(<p) which we know 
in terms of our initial data. Similarly we can determine -0-2 (4>) and ip-3{4>)- 
When we come to ^_ 4 (0) it seems that we would need the unknown ipo{<fi) 
which, fortunately, is not the case because ipo(4>) drops out of the evolution 
equation completely. It does not appear in the middle term on the left hand 
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side because now, for n = 0, VLi « — VLi /2_i = 0. Furthermore, we have seen 
as a general conclusion of loop quantizations that the matter Hamiltonian 
annihilates the singular state |0), which in the triad representation translates 
to -ff ma tter(" = 0) = independently of the kind of matter. Thus, -0o(0) drops 
out completely and ?/>_4 is determined solely by "04 • The further evolution to 
all negative n then proceeds without encountering any problems. 

Intuitively, we obtain a branch of the universe at times "before" the clas- 
sical singularity, which cannot be seen in the classical description nor in the 
Whcclcr-DeWitt quantization. Note, however, that the classical space-time 
picture and the notion of time resolves around the singularity; the system can 
only be described by quantum geometry. The branch at negative times col- 
lapses to small volume, eventually reaching volume zero in the Planck regime. 
There, however, the evolution does not stop, but the universe bounces to en- 
ter the branch at positive time we observe. During the bounce, the universe 
"turns its inside out" in the sense that the orientation of space, given by 
sgn(p), changes. 

In the discussion we ignored the fact that we could not determine ipo (4>) by 
using the evolution equation. Is it problematic that we do not know the values 
at the classical singularity? There is no problem at all because those values 
just decouple from values at non-zero n. Therefore, we can just choose them 
freely; they do not influence the behavior at positive volume. In particular, 
they cannot be determined in the above way from the initial data just because 
they arc completely independent. 

The decoupling of ipo(4>) was crucial in the way we evolved through the 
classical singularity. Had the values not decoupled completely, it would have 
been impossible to continue to all negative n. It could have happened that 
the lowest order coefficient is zero at some n, not allowing to determine ip n -4, 
but that this unknown value would not drop out when trying to determine 
lower ip n . In fact, this would have happened had we chosen a factor ordering 
different from the one implicitly assumed above. Thus, the requirement of a 
non-singular evolution selects the factor ordering in loop quantum cosmology 
which, in turn, fixes the factor ordering of the Wheeler DcWitt equation (13) 
via the continuum limit. One can then re-check results of Wheeler-DeWitt 
quantum cosmology which are sensitive to the ordering [12] with the one we 
obtain here. It is not one of the ordcrings usually used for aesthetic reasons 
such that adaptations can be expected. An initial step of the analysis has 
been done in [43]. 

To summarize, the evolution equation of loop quantum cosmology allows 
us, for the first time, to push the evolution through the classical singular- 
ity. The theory tells us what happens beyond the classical singularity which 
means that there is no singularity at all. We already know that energy den- 
sities do not diverge in a loop quantization, and now we have seen that the 
evolution does not stop. Thus, none of the conditions for a singularity is 
satisfied. 
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Dynamical initial conditions 

In the Wheeler-DeWitt quantization the singularity problem has been glossed 
over by imposing initial conditions at a = 0, which does have the advantage 
of selecting a unique state (up to norm) appropriate for the unique universe 
we observe. This issue appears now in a new light because n = does not 
correspond to a "beginning" so that it does not make sense to choose initial 
conditions there. Still, n = does play a special role, and in fact the behavior 
of the evolution equation at n = implies conditions for a wave function [46] . 
The dynamical law and the issue of initial conditions are intertwined with 
each other and not separate as usually in physics. One object, the constraint 
equation, both governs the evolution and provides initial conditions. Due to 
the intimate relation with the dynamical law, initial conditions derived in 
this way are called dynamical initial conditions. 

To see this we have to look again at the recurrence performed above. We 
noted that the constraint equation does not allow us to determine iio(4>) from 
the initial data. We then just ignored the equation for n = 4 and went on to 
determine the values for negative n. The n = 4-cquation, however, is part of 
the constraint equation and has to be fulfilled. Since -00 drops out, it is a linear 
condition between ipi and tps or, in a very implicit way, a linear condition for 
our initial data. If we only consider the gravitational part, i.e. the dependence 
on n, this is just what we need. Because the second order Wheeler-DeWitt 
equation is reproduced at large volume, we have a two-parameter freedom 
of choosing the initial values in such a way that the wave function oscillates 
only slowly at large volume. Then, one linear condition is enough to fix the 
wave function up to norm. When we also take into account the matter field, 
there is still more freedom since the dependence of the initial value on <f> is 
not restricted by our condition. But the freedom is still reduced from two 
functions to one. Since we have simply coupled the scalar straightforwardly 
to gravity, its initial conditions remain independent. Further restrictions can 
only be expected from a more universal description. Note also that there 
are solutions with a wave length the size of the Planck length which are 
unrestricted (since the evolution equation only relates the wave function at 
n and n ± 4). Their role is not understood so far, and progress can only be 
achieved after the measurement process or, in mathematical terms, the issue 
of the physical inner product is better understood. 

In its spirit, the dynamical initial conditions are very different from the 
old proposals since they do not amount to prescribing a value of the wave 
function at a = 0. Still, they can be compared at least at an approximate 
level concerning implications for a wave function. They are quite similar to 
DeWitt's original proposal that the wave function vanishes at a = 0. The 
value at a — itself would not be fixed, but quite generally the wave function 
has to approach zero when it reaches n = 0. In this sense, the dynamical 
initial conditions can be seen to provide a generalization of DeWitt's initial 
condition which does not lead to ill-posed initial value problems [47]. 
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For the closed model with k = 1 we can also compare the implications 
with those of the tunnelling and the no-boundary proposals which have been 
defined only there. It turns out that the dynamical initial conditions are very 
close to the no-boundary proposal while they differ from the tunnelling one 
[43]. 

5.4 Phenomenology 

So far, we have discussed mainly conceptional issues. Now that we know that 
loop quantum cosmology is able to provide a complete, non-singular descrip- 
tion of a quantum universe we can also ask whether there are observational 
consequences. We already touched this issue by comparing with the older 
boundary proposals which have been argued to have different implications 
for the likelihood of inflation from the initial inflaton values they imply. How- 
ever, the discussion has not come to a definite conclusion since the arguments 
rely on assumptions about Planck scale physics and also the interpretation 
of the wave function about which little is known. Loop quantum cosmology 
provides a more complete description and thus can add substantially to this 
discussion. An analysis analogous to the one in the Wheeler-DeWitt quan- 
tization has not been undertaken so far; instead a strategy has been used 
which works with an effective classical description implementing important 
quantum geometry effects and thus allows to sidestep interpretational prob- 
lems of the wave function. It provides a general technique to study quantum 
effects in a phenomenological way which is currently being used in a variety 
of models. 

Effective Friedmann equation 

The central idea is to isolate the most prominent effects of quantum geom- 
etry and transfer them into effective classical equations of motion, in the 
case of isotropic cosmology an effective Friedmann equation [48] . The most 
prominent effect we have seen is the cut-off observed for inverse powers of 
the scale factor (which can be thought of as a curvature cut-off). It is a non- 
perturbative effect and has the additional advantage that its reach can be 
extended into the semiclassical regime by choosing a large ambiguity param- 
eter j. 

In the equations for the isotropic model, inverse powers of the scale factor 
appear in the kinetic term of the matter Hamiltonian, e.g. (7) for a scalar. 
We have discussed in Section 2 that it is difficult to suppress this term by 
arranging the evolution of 4>. Now we know, however, that quantum geom- 
etry provides a different suppression mechanism in the inverse scale factor 
operator. This has already played an important role in showing the absence 
of singularities since in fact the matter Hamiltonian vanishes for n = 0. In- 
stead of a~ 3 we have to use a quantization of the inverse scale factor, e.g. 
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in the form a~ 3 whose eigenvalues (30) with m = 3 are bounded above. We 
can introduce the effect into the classical equations of motion by replacing 
d(a) = a -3 with the bounded function dj(a) = a~ 3 p(3a 2 / j^ip) 6 where p(q) 
is defined in (29) and we choose a half-integer value for j. The effective scalar 
energy density then can be parameterized as 

p eS {a) = \xa l ^-H- l(a) -*pl + W{cl>) 

with parameters x and I which also depend on quantization ambiguities, 
though not significantly. Note that £p appears in the denominator demon- 
strating that the effect is non-perturbative; it could not be obtained in a 
perturbative quantization in an expansion of G. With our function p{q) we 
have I = 12 for very small a (but it decreases with increasing a), and it is 
usually larger than 3, even taking into account different quantization choices, 
thanks to the high power of p(q) in dj(a). Thus, the effective equation of state 
parameter w = —1/3 in the parametrization (12) is smaller than — 1; quantum 
geometry predicts that the universe starts with an initial phase of inflation 
[48] . It is a particular realization of super- inflation, but since w increases with 
a, there is no pole as would be the case with a constant w. Note that this 
does not require any special arrangements of the fields and their potentials, 
not even an introduction of a special inflaton field: any matter Hamiltonian 
acquires the modified kinetic term such that even a vanishing potential im- 
plies inflation. Inflation appears as a natural part of cosmological models in 
loop quantum cosmology. Moreover, the inflationary phase ends automati- 
cally once the expanding scale factor reaches the value a « £p where 
the modified density reaches its peak and starts to decrease (Fig. 3). 




Fig. 3. Behavior of the scale factor (left) and the scalar field (right) during quantum 
geometry inflation (ending at t w 0.4 for j = 100) [43], both plotted in Planck units. 
The potential is just a mass term 8irGW(<p) = lO~ 3 ft0 2 /2, and initial conditions 
for the numerical integration are <f> a = 0, ^/K(j>o = lO -5 ^ 1 at ao = 2£p. 
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Inflation 

Thus, inflation appears as a natural consequence, but it is less clear what role 
it can play. For an inflationary period responsible for structure formation it 
has to last long enough (in terms of e-foldings, i.e. a large ratio of the final a 
and the initial a) and to be very close to standard inflation, i.e. w = — 1. T he 
final scale factor for quantum geometry inflation is easy to find, a w y/ j^j/Z lp 
as just discussed. The initial value, however, is more complicated. In a flat 
model, the inflationary period starts as close to a = as we want, but at those 
small values the effective classical description must break down. Moreover, in 
a closed model the region close to a = is classically forbidden 17 which also 
sets a lower limit for the initial a. All these issues depend more sensitively 
on the kind of matter added to the model and have not yet been analyzed 
systematically. 

An alternative application of quantum geometry inflation can be seen in 
combination with standard inflation. We have discussed that the standard 
scenario requires a special potential and also special, very large initial values 
for the inflaton. For instance, for chaotic inflation with the potential W{4>) = 
ir7j</> 2 we need to start with 4> > mp = h/lp which is huge compared to its 
own mass m. If we couple quantum geometry inflation with chaotic inflation, 
we would first observe an inflationary expansion at small volume which can 
stop at small a (i.e. j can be of the order one). During this phase also the 
evolution of the scalar is modified compared to the standard one since now 
dj(a) appears in the Hamiltonian equations of motion instead of a~ 3 . This 
leads to a differential equation 

^ dlogyq) ,._ a3 ^ (a) ^ w 

for <j>. For the always decreasing a~ 3 instead of dj(a) we obtain the previous 
equation (10) with the friction term. The modified dj(a), however, is increas- 
ing for small a such that we obtain a friction term with the opposite sign. 
This will require the inflaton to move up the potential, reaching large values 
even if it would start in (p(a = 0) = [43]; see Fig. 3. 

5.5 Homogeneous cosmology 

The framework of loop quantum cosmology is available for all homogeneous, 
but in general anisotropic models [35]. When we require that the metric of a 
homogeneous model is diagonal, the volume operator simplifies again allowing 
an explicit analysis [49] . One obtains a more complicated evolution equation 

17 There is, however, a mechanism which leads to a small classically allowed region 
for small a including a — even in a closed model [43]. This comes from a 
suppression of intrinsic curvature analogous to the cut-off of a -1 which would 
be a suppression of extrinsic curvature. 
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which is now a partial difference equation for three degrees of freedom, the 
three diagonal components of the metric. Nevertheless, the same mechanism 
for a removal of the classical singularity as in the isotropic case applies. 

This is in particular important since it suggests an absence of singulari- 
ties even in the full theory. It has been argued [50] that close to singularities 
points on a space-like slice decouple from each other such that the metric in 
each one is described by a particular homogeneous model, called Bianchi IX. 
If this is true and extends to the quantum theory, it would be enough to have 
a non-singular Bianchi IX model for singularity freedom of the full theory. 
Even though the classical evolution of the Bianchi IX model is very compli- 
cated and suspected to be chaotic [51], one can see that its loop quantization 
is singularity- free [52]. In fact, again the cut-off in the inverse scale factor 
leads to modified effective classical equations of motion which do not show 
the main indication for chaos. The Bianchi IX universe would still evolve in a 
complicated way, but its behavior simplifies once it reaches small volume. At 
this stage, a simple regular transition through the classical singularity occurs. 
This issue is currently being investigated in more detail. Also other homoge- 
neous models provide a rich class of different systems which can be studied 
in a phenomenological way including quantum geometry modifications. 

6 Quantum gravity phenomenology 

Since quantum gravity is usually assumed to hold at scales near the Planck 
length £ P ~ 10 -32 cm or, equivalently, Planck energy E P := h/£ P <~ 10 18 GeV 
experiments to probe such a regime were considered out of reach for most of 
the past. Recently, however, phenomena have been proposed which compen- 
sate for the tiny size of the Planck scale by a large number of small corrections 
adding up. These phenomena include in vacuo dispersion relations for gamma 
ray astrophysics [53, 54, 55], laser-interferometric limits on distance fluctua- 
tions [56, 57], neutrino oscillations [58], threshold shifts in ultra high energy 
cosmic ray physics [59, 60, 62, 61], CPT violation [63] and clock-comparison 
experiments in atomic physics [64]. They form the so called quantum gravity 
phenomenology [65]. 

The aim is to understand the imprint which the structure of space-time 
predicted by a specific theory of quantum gravity can have on matter prop- 
agation. Specifically, dispersion relations are expected to change due to a 
non-trivial microscopic structure (as in condensed matter physics where the 
dispersion relations deviate from the continuum approximation once atomic 
scales are reached). For particles with energy E <C Ep and momentum p the 
following modified vacuum dispersion relations have been proposed [54] : 

p 2 = E 2 (l+£E/E P +0((E/E P ) 2 )) , (33) 

where £ ~ 1 has been assumed, which still has to be verified in concrete 
realizations. Furthermore, this formula is based on a power series expansion 
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which rests on the assumption that the momentum is analytic at E = as 
a function of the energy. In general, the leading corrections can behave as 
(E/E P ) r+1 , where T > is a positive real number. 

Since the Planck energy is so large compared to that of particles which 
can be observed from Earth, the correction would be very tiny even if it is 
only of linear order. However, if a particle with the modified dispersion rela- 
tion travels a long distance, the effects can become noticeable. For instance, 
while all photons travelling at the speed of light in Minkowski space would 
arrive at the same time if they had been emitted in a brief burst, Eq. (33) 
implies an energy dependent speed for particles with the modified dispersion 
relations. Compared to a photon travelling a distance L in Minkowski space, 
the retardation time is 

Atm £LE/Ep . (34) 

If L is of a cosmological scale, the smallness of E/E P can be compensated, 
thus bringing At close to possible observations. Candidates for suitable sig- 
nals are Gamma Ray Bursts (GRB's), intense short bursts of energy around 
E ~ 0.20MeV that travel a cosmological distance L ~ 10 10 ly until they reach 
Earth. These values give At ~ 0.01ms which is only two orders of magnitude 
below the sensitivity St for current observations of GRB's [66, 67] (for planned 
improvements see [68] ) . For the delay of two photons detected with an energy 
difference AE, the observational bound £p/£ > 4x 10 16 GeV was established 
in [69] by identifying events having AE = 1 TeV arriving to Earth from the 
active galaxy Markarian 421 within the time resolution At = 280 s of the 
measurement. Moreover, GRB's also seem to generate Neutrino Bursts (NB) 
in the range 10 5 — 10 10 GeV in the so-called fireball model [70, 71] which can 
be used for additional observations [72, 58, 73]. 

In summary, astrophysical observations of photons, neutrinos and also 
cosmic rays could make tests of quantum gravity effects possible, or at least 
restrict possible parameters in quantum gravity theories. 

Within loop quantum gravity attention has focused on light [74, 75] and 
neutrino propagation [76]. Other approaches aimed at investigating similar 
quantum gravity effects include string theory [77] , an open system approach 
[78], perturbative quantum gravity [79, 80] and non commutative geome- 
try [81]. A common feature to all these approaches is that correction terms 
arise which break Lorentz symmetry. These studies overlap with a system- 
atic analysis providing a general power counting renormalizable extension of 
the standard model that incorporates both Lorentz and CPT violations [82] . 
Progress in setting bounds to such symmetry violation has been reported in 
[55, 64, 83, 84, 85]. 

6.1 An implementation in loop quantum gravity 

In order to implement the central idea, one needs states approximating a 
classical geometry at lengths much larger than the Planck length. The first 
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proposed states of this type in loop quantum gravity were weave states [86] . 
Flat weave states \W) with characteristic length C were constructed as in 
Section 4 using collections of circles of Planck size radius (measured with the 
classical background geometry to be approximated) in random orientation. 
At distances > £ the continuous flat classical geometry is reproduced, 
while for distances d <C C the discrete structure of space is manifest. The 
search for more realistic coherent states, which not only approximate the 
classical metric, but also its conjugate, the extrinsic curvature, is still ongoing 
[87, 88, 89, 90]. Current calculations have been done at a heuristic level by 
assuming simple properties of semiclassical states. The prize to pay is that 
some parameters, and even their order of magnitude, remain undetermined. 
Thus, these studies explore possible quantum gravity effects without using 
details of specific semiclassical states. Once properties of semiclassical states 
become better known, one can then check if the existing calculations have to 
be modified. 

The setup requires to consider semiclassical states \S) for both gravity, 
the background for the propagation, and the propagating matter. One has 
to require that they are peaked at the classical configurations of interest 
with well defined expectation values such that there exists a coarse-grained 
expansion involving ratios of the relevant scales of the problem. Those are 
the Planck length £p, the characteristic length C and the matter wavelength 
A satisfying £p <C C < A; see Fig. 4. The effective Hamiltonian is thus defined 

by 

^Matter : = < & I -^Matter | S > . (35) 




Fig. 4. Representation of the different scales of the problem. A coarse graining 
scale C ~ A indicates Planck length features £ P are minute as compared to matter 
scales. L here represents the size of a large piece of space. 
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Light 

The full quantum Hamiltonian for the electromagnetic field is of the form 
(20). With the above assumptions about semiclassical states one can arrive 
at an effective electromagnetic Hamiltonian [75] 



B 2 + E 2 



i(i + 7 (£ P /cf +2T ) ( : 
+ e 3 £l ( b ■ V 2 B + E ■ V 2 E) 

+6 8 l P (B_- (V x B) + E- (V x E)) +■■■], (36) 

up to order £\ and neglecting non linear terms. The coefficients Qi have not 
yet been derived systematically; rather, the expression is to be understood 
as a collection of all the terms which can be expected. Precise values, and 
even the order of magnitude, can depend significantly on the explicit proce- 
dure followed to obtain the values from a semiclassical state. Moreover, as 
usual there are quantization ambiguities in the quantum Hamiltonian which 
influence the coefficients of the correction terms [41]. 

From the effective Hamiltonian (36) we obtain the equations of motion 

BE 

A(V x B) — ■= + 2£l6 3 V 2 (V x B) - 29 S £ P V 2 B = 0, (37) 
A(V xE) + -^ + 24<9 3 V 2 (V x E) - 26 8 e P V 2 E = 0, (38) 



where 



1 + 7 (£ P /C) 2+2T . (39) 



The above equations are supplemented by V • B_ = 0, together with the 
constraint V • E_ = 0, appropriate for vacuum. 

Modifications in the Maxwell equations (37) and (38) imply a modified 
dispersion relation which, neglecting the non-linear part, can be derived by 
introducing the plane wave ansatz 

E = E c l{kx - Ut \ B = B c l{kx - Ut \ k=\k\. (40) 

The result is 

uj = k (l + 7 (£ P /C) 2+2T - 2 6 3 (Hp) 2 ± 26> 8 (H P )) (41) 

where the two signs of the last term correspond to the different polarizations 
of the photon. The speed of a photon becomes 



dk 



fc = 1 ± 4 6 8 (k£ P ) - 66 3 (k£ P ) 2 + 6 7 {k£ P ) 2+zl + ■■■ (42) 
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The scale C has been estimated by its maximal value 1/k. Clearly (42) is 
valid only for momenta satisfying (^p fc) « 1. 

There are also possible non-linear terms in the effective Maxwell equations 
[75]. They can become significant in strong magnetic fields, but the correc- 
tions obtained in the corresponding refraction indices are much smaller than 
similar effects in Quantum Electrodynamics. Nevertheless, quantum gravity 
corrections have distinct signatures: a main difference is that the speed of 
photons with polarization parallel to the plane formed by the background 
magnetic field and the direction of the wave is isotropic. 



Spin- 1/2 particles 

Similarly, one can derive an effective Hamiltonian for a spin-i field of mass 
m [76]: 



spin 



+ 



i J d 3 x [ir(x)T d d d A £(x) + c.c. + {AC)- 1 n(x)C£(x) (43) 
a 2 (a + f3tpr a d a ) £(x) + tt t (x) (a + (3i P T a d a ) <t 2 tt(x)}}, 



where 

A = (i + K! (e P /cf +1 + k 2 (£ P /cf r+2 + ± k 3 4 v 2 ) 

C = \k 7 (l P /£) r 1% V 2 

a = 1 + k 8 (£p/C) t+1 , /3 = k 9 + Kn (e P /C) r+1 
This leads to wave equations 



(44) 



d i „ . C 
— - A a- ■ V - i— 
dt 2C 



9 i „ . C 
dt +A(T - V + l 2C 



£(t, x) + m(a- <x • V) x (t, x) = (45) 

X (t, x)+m(a- \i(3£ v ct • V) f (i, x) = (46) 



with x(t, x) — i a 2 C* (t, x). As before, the dispersion relation can be obtained 
by inserting plane wave solutions, this time positive and negative energy 
solutions 

W(p, h)e^ iEt/R±ip - x/h (47) 
where W(p,h) are helicity (ct ■ p) eigenstates, with h = ±1, so that 



-e' 1 ' 1 ' sin 



W 2 A (4*) 
cos(0/2) J ■ [ 8J 



For ultra-relativistic neutrinos (p » to) one obtains 
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£ P E ± {p, C) = p£ P + £ P m 2 /2p ± \ (£ P m) 2 k 9 - \ k 3 (£ P pf 
+ (£ P /C) r+1 \k iP £ p T jK7 (£pp) 2 } + {tp/C) 



2T+2 



K2P£p 



(49) 



and 




v± (p, C) = 1 - ^ - § S3 (^PP) 2 + ( V£) T+1 («i T § Mpp) + {£p/C) 2r+2 K 2 . 



There are two physically interesting effects related to the dispersion rela- 
tions just described for neutrinos. Namely neutrino oscillations for different 
flavors and time delay between neutrinos and photons coming from the same 
GRB. Estimates in this respect have been obtained in [76]. 

6.2 Summary 

The phenomenological considerations described here are intended to give an 
idea of possible consequences of quantum gravity corrections. They start with 
assumptions about a state approximating a classical flat metric, a classical 
flat gravitational connection and a generic classical matter field, at scales 
larger than the coarse-grained characteristic length C » £p. Under these 
assumptions, modified dispersion relations can be expanded in the Planck 
length. In general, there are different types of corrections, which can have 
different dependence on, e.g., the helicity or the scale C. This also includes 
the parameter Y encoding our (current) ignorance of the scaling of the grav- 
itational connection in a semiclassical expectation value. 

The following motivations for the value of Y have been made [75]: (i) 
Y = can be understood as that the connection can not be probed below 
the coarse graining scale C. The corresponding correction scales as (k£p) 2 . 
(ii) Y = 1 may be interpreted as the analog of a simple analysis [91], based 
on a saturation of Heiscnberg's uncertainty relation inside a box of volume 
C 3 : AE - £ P /C, A A - £ P /C 2 and AEAA - Gh/C 3 . Then the correction 
behaves as (k£p) 4 . (iii) A value Y = — | would lead to a helicity independent 
first order correction (i.e. (k£p))\ a negative value, however, is not allowed. 
Further fractional values have been obtained in [92] from a detailed proposal 
for coherent states in loop quantum gravity [87, 88]. 

From an observational point of view, lower order correction terms would 
certainly be preferable. Most of the evaluations so far have been done for first 
order terms, but recently also higher order corrections have been started to 
be compared with observations [62, 93]. 



7 Outlook 



As discussed in this article, loop quantum gravity is at a stage where physical 
results are beginning to emerge which will eventually be confronted by obser- 
vations. To obtain these results, as usually, approximation schemes have to be 
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employed which capture the physically significant contributions of a full the- 
ory. In our applications we used the minisuperspace approximation to study 
cosmological models and a semiclassical approximation for the propagation 
of particles. We have to stress, however, that these approximation schemes 
are currently realized at different levels of precision, both having open issues 
to be filled in. Loop quantum cosmological models are based on symmet- 
ric states which have been explicitly constructed as distributional states in 
the full theory. There are no further assumptions besides the central one of 
symmetries. A partially open issue is the relation of symmetric operators to 
those of the full theory. A precise derivation of this relation will complete our 
understanding of the models and also of the full theory, but it is not expected 
to imply changes of the physical results since we know that they are robust 
under quantization ambiguities. 

As for loop quantum gravity phenomenology its central ingredient are 
semiclassical states which are being investigated with different strategies 
leading to different proposals. The present explicit calculations are based 
on simple assumptions about semiclassical states which have to be probed in 
an eventual realization. Thus, there is not only a central simplifying assump- 
tion, semiclassicality, but also additional assumptions about its realization. 
These assumptions affect the presence as well as the magnitude of possible 
correction terms. It is not just the relation between phenomenology and ob- 
servations, but also the one between phenomenology and the basic theory 
which has to be understood better. 

Furthermore, there are important conceptual issues which are not yet 
completely understood. For instance, it would be essential to see the emer- 
gence of a classical space-time from semiclassical quantum states in order to 
study a particle moving in a state which approximates Minkowski space. A re- 
lated issue is the fact that the discreteness of quantum geometry is supposed 
to lead to correction terms violating Lorentz symmetry. Such a violation, in 
turn, implies the existence of a distinguished time- like vector. An open con- 
ceptual issue is how such a distinguished vector can arise from the discrete 
formulation. 18 For this purpose one would need a distinguished rest frame 
which could be identified using the cosmic background radiation [74]. 

Future work will progress along several lines according to the different 
open problems. First, at a basic level, the conceptual issues will have to be 
understood better. In the case of quantum gravity phenomenology this will 
come as a consequence of additional insights into semiclassical states which 
are under investigation [87, 88, 89, 90]. This will also change the way how 
explicit calculations are implemented, and the precision of known results will 
be enhanced leading to a stronger confrontation with observations. Finally, 
there are many phenomenological effects which have not yet been investigated 

18 Models to understand modifications of the usual Lorentz symmetry have been 
developed in [94, 95]. 



Loop quantum gravity 39 



in the context of loop quantum gravity. Loop effects will lead to changes whose 
significance regarding observations has to be studied. 

Already the present stage of developments proves that loop quantum grav- 
ity is a viable description of aspects of the real world. It offers natural solu- 
tions to problems, as e.g. the singularity problem, which in some cases have 
been open for decades and plagued all other theories developed so far. At 
the same time, sometimes surprising consequences emerged which lead to a 
coherent picture of a universe described by a discrete geometry. All this es- 
tablishes the viability of loop quantum gravity, and we are beginning to test 
the theory also observationally. 
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